Dynamical dipole 7 radiation in heavy-ion collisions on the basis of the quantum 

molecular dynamics model 
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Dynamical dipole 7-ray emission in heavy-ion collisions is explored in the framework of the quan- 
tum molecular dynamics model. The studies are focused on systems of ''"Ca bombarding ^*Ca and 
its isotopes at different incident energies and impact parameters. Yields of 7 rays are calculated 
and the centroid energy and dynamical dipole emission width of the 7 spectra are extracted to in- 
vestigate the properties of 7 emission. In addition, sensitivities of dynamical dipole 7-ray emission 
to the isospin and the symmetry energy coefficient of the equation of state are studied. The results 
show that detailed study of dynamical dipole 7 radiation can provide information on the equation 
of state and the symmetry energy around the normal nuclear density. 

PACS numbers: 25.70.Gh, 24.30.Cz, 25.70.Ef, 25.70.Lm 



Giant dipole resonance (GDR) 7-ray emission built on 
heavy-ion collisions, which was predicted in early macro- 
scopic models (T], Q has been widely studied both exper- 
imentally and theoretically during the past few decades 
d, 0) ■ It can be ascribed to a collective vibration of pro- 
tons against neutrons with a dipole spatial pattern that 
starts with the overlap of two nuclei in both fusion or in- 
complete fusion. Recently, a so-called dynamical dipole 
mode (pre-equilibrium GDR) that is generated by a high- 
amplitude dipole collective motion of nucleons before the 
formation of a fully equilibrated compound nucleus (CN) 
has been discussed [5|-l7| and it can encode information 
about the early stage of collisions when the CN is still 
in a highly deformed configuration. Some efforts have 
been made to study the dipole resonance formed during 
fusion in N /Z asymmetry heavy- ion reactions as it could 
be a good probe for the symmetry term in the equation 
of state (EOS), which is one of the hotly debated issues 
related to nuclear astrophysics problems such as the ele- 
ments burning in supernova Q. 

Several microscopic transport simulations, such as 
time-dependent Hartree-Fock [H, Boltzmann-Nordheim- 
Vlasov [9!] , and constrained molecular dynamics , have 
been successfully used to study GDR dynamical dipole 
emission. However, the quantum molecular dynamics 
(QMD) model, which has been effectively used to dispose 
problems of heavy- ion collisions up to 2 GeV/nucleon 
has not been used to investigate GDR. In princi- 
ple, energy is not conserved in the QMD model and the 
model is more suited for higher-energy phenomena. But, 
it is still of interest to see how the QMD model works for 
GDR 7 radiation and discuss its properties . 

We use the isospin-dependent QMD (IDQMD) model, 
which is based on the general QMD model but embod- 
ies the isospin degrees of freedom in two-body nucleon- 
nucleon collisions and Pauli blocking, as well as consid- 



ering the difference between neutron and proton density 
distributions for nuclei far from the /3-stability line and 
treating the sampling of neutrons and protons in the ini- 
tialization in phase space, respectively [T2|. [isj. 



The mean field involved in the present IDQMD model 
is given by U{p) ^ U^^^ + C/^°"' + C/'^"'^ + W^"", with 
[/Sky^ jjConX^ jjYuk ^^^^ jjsym ^presenting the Skyrme 

potential, the Coulomb potential, the Yukawa potential 
and the symmetry potential interaction, respectively 
The Skyrme potential is U^^^ = a{p/po) + i3{p/poy, 
where po = 0.16 fm~^(the saturation nuclear density) 
and p is the nuclear density. We take the parameters 
a = -356 MeV, /3 = 303 MeV, and 7 = 7/6, correspond- 
ing to a soft EOS, and a = -124 MeV, /3 = 70.5 MeV, 
and 7 = 2, corresponding to a stiff EOS. JJ^^^ is a long- 
range interaction that is necessary to describe the surface 
of the nucleus and takes the following form: fj^"*^ — 
(yy/^)Y^i^j ^xpiLm'^)/rij ■ [exp{mrij)erfc{VLm - 

Tijjsf^V) — exp{mrij)erfc{\/Lm + rij/\/TL)] with 



Vy = 0.0024GeV, 



= 0.83, and L is the so-called 
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Gaussian wave-packet width (here L = 2.0 fm ). r,j is 
the relative distance. In principle, if one adds JJ^'^^ into 
the effective potential, some part of the two-body term 
in C/'°'^ should be subtracted. Therefore a should be 
smaller [11]. However, considering that we aim to discuss 
some qualitative effects of the stiffness of potential on 
pre-equilibrium GDR 7 emission, we tentatively keep 
the preceding two sets of potential parameters. 



In recent years, some efforts have been made to con- 
strain the symmetry energy term of the EOS from nuclear 
reactions induced by asymmetric heavy ions, and vari- 
ous probes have been proposed Q . Considering that the 
pre-equilibrium GDR 7 ray originates from the isospin 
nonequilibrium of the entrance channel or the composite 
intermediate system, its spectra should be dependent on 
the symmetry potential. The symmetry potential in the 
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IDQMD model is obtained by 

C*gvm X ^ 1 



jjsym 



2po 



(47rL)3/2 



exp[- 



4L 



(1) 



with Cgym the symmetry energy strength, being the 
zth component of the isospin degree of freedom, which 
equals 1 or —1 for neutrons or protons, respectively. 
A detailed description of the QMD model is given in 
Ref.[ll|. 

The giant dipole moment in coordinator space [DR{t)] 
and in momentum space [DK(t)] is written, respectively, 
as follows [145 : 

NZ NZ 
DR{t) ^ —X{t) ^ —{Rz - Rn), (2) 

where X{t) is the distance between the centers of mass of 
protons and neutrons, Rz = ^iip) ^^'^ Rn =Y1 ^ii^^) 

i i 

is the center of mass of protons and neutrons, respec- 
tively, and 

, , NZ T-r, ^ NZ P„, 



is just the canonically conjugate momentum of the 
DR{t), where J^(t) denotes the relative momentum, with 
Pn (Pp) centers of mass in momentum space for neutrons 
(protons). 

Here N = Np + Nt a.nd Z ^ Zp + Zt, the suffix P 
(or T) marks projectile (or target). We have adopted the 
initial dipole moment as follows: 



NZ 

DR{t = 0) = —\Rz{t = 0) - RN{t = 0)1 
= ^ ^ — ^ZpZt\{ — )t - (^)-pI' 



(4) 



where Rp and Rt are the radii of the projectile and 
target, respectively. 

Derived from the overall dipole moment I?(i), we can 
get the 7-ray emission probability for energy E Q 



dP 
dE 



,dVk 



(5) 



Stt Ehc dt 

where dP/dE can be interpreted as the average number 

of 7 rays emitted per energy unit, while ^^{E) is the 
Fourier transformation of the time derivative on the kth 
components x and z of 



(6) 



where V denotes the dipole resonance of the system and 
was derived from the total dipole according to the equa- 
tion 



jr dD{t)k 



(7) 
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FIG. 1: (Color online) Left: Evolution of the giant dipole mo- 
ment with time in coordinate space [Dx{t)] and momentum space 
[DKx{t)]. Right: Dynamical dipole mode [Vi(t)]. The reaction 
system is *'^Ca-|-^*Ca collision at 10 MeV/nucleon with 6= 1 fm 
and soft EOS parameters. 
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FIG. 2: (Color online) Left: Experimental data (stars) llOll to- 
gether with our IDQMD model [filled (red) circles] of the 7-ray 
spectrum of the ^''Ca-l-^^Ca system at 10 MeV/nucleon. Right: 
Measured value of 7-ray yield for the ^^0-|-^^®Sn system at 15.6 
MeV/nucleon (open circles) and our calculation results [filled 
(blue) triangles] under the same conditions. 



Figure [T] shows the evolution of the giant dipole mo- 
ment with time in coordinate space and momentum space 
for 4"Ca + ^^Ca colhsions at 10 IVIeV/nucleon and 6=1 
fm. The giant dipole moment attenuates rapidly after 
t — 400 fm/c, illustrating a clear reduction of the col- 
lective behavior. The similar behavior is observed for 
the pre-equilibrium or dynamical dipole mode [denoted 

To determine the reliability of our calculations of 
dipole emission with the IDQMD model, experimental 
data for the '^^Ca+^^Ca system [13] at £;beam = 10 
MeV/nucleon and the ^'^O-h^^Sn system 0| at i^bcam = 
15.6 MeV/nucleon are compared with our results in Fig. 
[21 It should be noted that the experimental data on 
dynamical dipole emission here are obtained by taking 
away the statistical contribution (results of the statistical 
model) from the measured spectra. In addition, calcula- 
tions for both systems are carried out with soft Skyrme 
potential and impact parameters being a triangle distri- 
bution between 1 and 6 fm. The good agreement between 
the data and the IDQMD calculation encourages us to 
make systematic calculations on pre-equilibrium GDR 7 
emission, such as its dependences on impact parameter, 
isospin, EOS, etc. 

So far both experimental data and theoretical calcula- 
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FIG. 3: (Color online) Left: Dipole dynamical emission spectra 
of the -^^Ca+^^Ca system at 10, 15, and 20 MeV/nucleon, with 
impact parameters 6 = 1 fm and soft EOS parameters. Right: 
Dynamical dipole 7-ray yield distribution of the ^''Ca+^^Ca system 
at 10 MeV/nucleon in the early 400 fm/c, with impact parameter 
b varied from 1 to 5 fm and soft EOS parameters. 
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FIG. 4: (Color online) Left: Central energy Ec of the system 
*''Ca+^*Ca and corresponding maximum 7-emission probability 
{dP/dE-y)x as a function of b with i?beam = 10 MeV/nucleon and 
soft EOS parameters. Right: Corresponding variation of the pre- 
equilibrium GDR, width of the system. 



tions have shown that both the centroid and the width 
of the pre-equihbrium dipole component, with the excep- 
tion of the 7-ray yield of GDR, are almost independent 
of the beam energy, and in this study, a similar con- 
clusion was drawn from the left panel in Fig. [31 where 
we report calculations of the dynamical 7-ray emission 
of ^oCa-h-^^Ca at -Bbeam = 10, 15, and 20 MeV/nucleon, 
respectively. In the figure, the centroid and width of 7 
spectra seem insensitive to the incident energy, but the 
7-ray emission probability dP/dE is suppressed at higher 
incident energies. A direct reduction of the initial charge 
asymmetry owing to pre-equilibrium nucleon emission 
(mostly neutrons) may account for the depressed 7 yield. 
In addition, a higher excitation energy deposited in the 
composite system can damp the GDR emission. The sup- 
pressed behavior of 7-ray emission with increasing inci- 
dent energy means that the heavy-ion collisions at a rel- 
ative low energy are suitable for studying the collective 
dipole resonance. In fact, the conclusion here was sup- 
ported by the existence of the dynamical dipole mode 
studied in fusion [T5l - [T7| and deep-inelastic heavy-ion col- 
lisions [Ulil. 

The right panel in Fig. [Hlgives the 7-ray emission prob- 
ability of the '**'Ca-|-'*^Ca system at 10 MeV/nucleon on 
a scale of 0-400 fm/c, with h = 1-5 fm and soft EOS 
parameters. Both the yield and the central energy of 
the 7-ray distribution show a decreasing tendency with 
a larger impact parameter. As we know, a larger im- 
pact parameter inevitably induces a larger time scale of 
the isospin variation, corresponding to a low frequency 
of the dipole resonance, and thus leads to a lower central 
energy. Besides, collisions with different impact param- 
eters, corresponding to different angular momenta, may 
generate dinuclei of different geometrical shapes [1^ , and 
necessarily, a larger impact parameter produces a dinu- 
clei system with a larger deformation. Using the conclu- 
sion in Ref. JJi] , the GDR in a deformable CN is expected 
to have a central energy lower than that in a spherical 
nucleus of similar mass. Therefore the suppression of the 
central energy with an increasing impact parameter can 
also be partly attributed to the larger deformation of the 
dinuclei at the emission moment. What is more, the de- 



formation of the CN may also bring about an anisotropic 
angular distribution pattern of the 7 ray, which implies 
that investigation of the anisotropic angular distribution 
of GDR emission may also be a promising way to obtain 
information on the deformation of the composite system, 
which deserves further investigation. 

Dependences of the centroids of pre-equilibrium GDR 
7-ray spectra and maximum yields on the impact pa- 
rameter are plotted in the left panel in Fig. 4; both 
the central energy (i.e., peak energy, denoted Ec) and 
the corresponding maximum 7-ray yield, {dP/dEj)x, de- 
crease monotonously with increasing impact parameter. 
Meanwhile, the slower isospin equilibration process for a 
large impact parameter may enhance the emission of low- 
energy and multicomponent 7 rays, and this could ex- 
plain why the dynamical dipole emission width (denoted 
FWHM) increases with the impact parameter, which is 
shown in the right panel in Fig|31 

The sensitivity of the 7-ray yield to the stiffness of the 
EOS is also studied by comparing the simulation with 
soft and stiff Skyrme potentials. The evolution of a giant 
dipole moment for the ^°Ca-|-^^Ca system, at i?bcam = 10 
MeV/nucleon and 6 = 4 fm, with different EOS parame- 
ters, soft EOS (filled circles) and stiff EOS (open stars), 
are plotted in the left panel in Fig\5\ the corresponding 7- 
emission probability of the system, in the right panel. As 
shown in the figure, the stiff Skyrme potential generally 
leads to larger centroid and stronger dynamical dipole 7 
ray than the soft one. 

Furthermore, the correlation between the dynamical 
dipole emission and the symmetry term of the EOS is 
also investigated by changing the symmetry energy coef- 
ficient (Csym)- In the left panel in Figl6l the spectrum of 
dynamical emission of the same system is reported, and 
the calculation is performed in the same situation except 
for the Csym- We can clearly see an increased yield of 7 
rays for larger Csym values. 

As mentioned before, another interesting aspect of the 
GDR emission property is its dependence on the N/Z 
asymmetry of the reaction systems, which in fact is a 
direct result of Eq.(6). This relationship, indeed, is of 
most importance for its association with the symme- 
try term of the EOS. Presently, taking the charge and 
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FIG. 5: (Color online) Left: The giant dipole moment evolves with 
time for the *''Ca+''*Ca system at i?beam = 10 MeV/nucleon and 
6 = 4 fm; Right: The corresponding 7-emission probability of the 
system. Symbols are defined in the keys. 
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FIG. 6: (Color online) Left: Symmetry energy coefficient (Csym) 
dependence of 7-emission probability for the system '"'Ca+''*Ca at 
^bcam = 10 MeV/nucleon and b = 4fm. Right: 7-emission prob- 
ability of systems '*'^Ca+^-^Ca, ^''Ca+'**Ca and reference system 
""Ca-t-^OCa at Ebcam = 10 MeV/nucleon and 6 = 4 f m with soft 
EOS (Csym = 32 MeV). Inset: Ratio of 7 yields of systems with a 
*^Ca target or ^-^Ca target to those with a '"^Ca target. Symbols 
are defined in the keys. 



mass symmetric system '*'^Ca-|-^°Ca as a reference sys- 
tem, whose pre-equilibrium effect is not expected because 
of small initial dipole moment, we performed calculations 
of dipole 7 emission in N /Z asymmetry systems, namely, 
^'^Ca-H^'^Ca and ^^C&+'^^Ca.. In the right panel in Fig. 
6, we show our results for 7 emission in the early 400 
fm/c at £'bcam — 10 MeV/nucleon, with = 4 fm and 
soft EOS. In the inset, the ratios of 7 yields are also 
shown, where the ratio between the '*°Ca-f'*^Ca and the 
''°Ca-|-^°Ca systems is shown by filled circles and the ratio 
between the ^°Ga.+^'^Ca and the **°Ca -I- *^°Ca systems is 



denoted by open triangles. From the figure, we find that 
the N/Z asymmetry systems ''°Ca-|-52cg^ (^/^ ^ i ^^-^ 

and *°Ca+"'*Ca {N/Z — 1.2) show a clear enhanced yield 
of the order of about 35% and 20%, with respect to the 
symmetry system '^°Ca-|-^°Ca {N/Z = 1.0), consistent 
with results in Ref.[3| and references therein. 

Experimental observations in Refs.fl^and [13] demon- 
strated that statistical 7 emission is insensitive to a par- 
ticular reaction entrance channel when we select the same 
phase-space region of the initial CN. Therefore, the evi- 
denced extra 7 yield can be explained by the dynamics 
of the compound system formed by a target and pro- 
jectile with a large N/Z asymmetry. In contrast, the 
equilibrium GDR width and the peak energy are rather 
insensitive to isospin asymmetry 2li]. 

In summary, we have applied the QMD model to inves- 
tigate pre-equilibrium giant dipole oscillations (the dy- 
namical dipole) and made systematical calculations for 
the low-energy reactions of ^°Ca+Ca isotopes. The re- 
sults show that dynamical dipole 7 radiation is a good 
probe for gathering important information on the early 
stage of fusion, in particular, the charge asymmetry en- 
trance channel. The dependences of the beam energy and 
impact parameter of the dipole resonance are discussed. 
The extra yield of pre-equilibrium GDR emission for sys- 
tems with with a large N/Z asymmetry suggests a strong 
dependence of dynamical emission on isospin asymmetry 
of the CN or, in other words, on the symmetry term of 
the EOS of nuclear matter. Meanwhile, the sensitivi- 
ties of dynamical spectrum properties to the mean field 
potential and its symmetry energy portion indicate that 
detailed study of dynamical dipole 7 radiation can be 
used as a probe to study the EOS around the normal 
nuclear density. 
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